Integrable multi atom matter-radiation models without rotating wave approximation 
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Interacting matter-radiation models close to physical systems are proposed, which without ro- 
tating wave approximation and with matter-matter interactions are Bethe ansatz solvable. This 
integrable system is constructed from the elliptic Gaudin model at high spin limit, where radiative 
excitation can be included perturbatively. 
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In dealing with interacting matter-radiation systems 
the counter rotating wave (CRW) terms, which inevitably 
appear in general situation, are usually neglected in- 
voking the rotating wave (RW) approximation, with- 
out which physical models generally become unsolvable. 
However the RW approximation breaks down away from 
the resonance condition at w„ k as well as for high 
intensity fields. Consequently, fast oscillations (with fre- 
quency uj a + ujf) can no longer be neglected compared to 
the slow ones (with frequency \u> a — u)f\) and one is forced 
to consider the general case with CRW terms. Moreover, 
since the situation is generic, this problem can arise in 
diverse models like those in quantum optics, in cavity 
QED both in microwave and optical domain [1, 2], in 
trapped ions irradiated by laser beams [3] as well as in 
transport through quantum dots coupled to boson model 
[4] and in quantum information transfer protocol with a 
superconducting circuit [5]. 

Under RW approximation one can obtain exactly solv- 
able models like Jayenes-Cummings (JC) [6], Buck- 
Sukumar [7] model and their multi-atom generaliza- 
tions [8-10] as well as g-deformed matter-radiation mod- 
els, inducing anisotropic together with higher nonlin- 
earity [10, 11]. However when this approximation is 
avoided, CRW terms appear having the form like H crw = 
P(b^a + + ba~) , in the simplest case and generally spoil 
the solvablity of the system. Models with CRW terms in 
various forms were investigated earlier [3, 12, 13], though 
such exactly solvable multi-atom models which are close 
to physical systems are not known in the literature, ex- 
cept perhaps a recent proposal [14] involving some un- 
physical terms. 

We propose here integrable multi-atom JC type 
matter-radiation models which include CRW terms and 
allow exact Bethe ansatz solutions for the vanishing radi- 
ation frequency, when the RW approximation is naturally 
not applicable. Field excitation term can be considered 
over the integrable model, perturbatively taking w/ to 
be small, while the interatomic interactions through spin- 
spin coupling can be included in an exact way. We derive 



our model from the elliptic Gaudin model [15], through 
spin-i representation for the atoms and the bosonic re- 
alization under high spin limit for the single mode ra- 
diation field. For following the logic of the construc- 
tion, we recall that all known exactly solvable multi- 
atom matter-radiation models with RW approximation 
are derived from integrable xxx, xxz spin models [8- 
10] or from their corresponding Gaudin limits [16, 17]. 
Note that at higher spins the underlying algebras associ- 
ated with these models are either su(2), su(l, 1) or their 
quantum deformations, both allowing the crucial bosonic 
realization. For possible construction of integrable mod- 
els without RW approximation, one may therefore expect 
to repeat the same construction starting from a more 
general integrable inhomogeneous xyz model with higher 
spin representation [18, 19]. The representative Lax op- 
erator of this model depends, apart from the spectral 
parameter A and the anisotropy parameter a, on the el- 
liptic modulus k, inhomogeneity parameters z n as well 
as on the spin s representation, through its dependence 



on coefficients W P (X 



y, k), p = 0, 1, 2, 3, expressed 



through elliptic functions and operators S p (s,a,k) sat- 
isfying the Sklyanin algebra [20]. In fact, all known in- 
tegrable matter-radiation models mentioned above can 
be obtained from this general structure at different lim- 
its of the parameters involved. For example, at k = 0, 
one recovers the trigonometric xxz case, when generators 
of the Sklyanin algebra reduce to those of the quantum 
algebra: S a (s,a),a — 1,2,3. Under a further limit of 
a — > one obtains the trigonometric Gaudin model with 
the operators reducing to S a (s), satisfying the standard 
Lie algebra. If however at the same time a limit A — » is 
imposed on the spectral parameter (together with z n ), we 
recover the XXX C0jSG 0jS well as the corresponding rational 
Gaudin model. 

We however keep the elliptic modulus k nontrivial 
throughout our construction and thereby follow a rout 
different from the earlier ones. Observe that, unlike 
the above cases the Sklyanin algebra does not have 
any known bosonic realization, though fortunately at 
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a — > it reduces to the standard algebra allowing bosonic 
mapping through the Holstein-Primakoff transformation 

(HPT): S+ = tfy/s-pN, S~ = y/s - pNb, S 3 = 
N — p§, N = b^b, where p = ± correspond to su(2) 
(su(l,l)). At this limit the xyz model reduces to the 
integrable elliptic Gaudin model [15] with mutually com- 
muting conserved quantities 

H n = Y, ™a{nm)S a n S? nl (1) 

with Wl (X) = f|(A), w 2 (X) = f£(A), w 3 (X) = £(A) 
and H n = 0, where we have introduced short hand 
notation w a (nm) = w a (z n — z m ). However, we face now 
the difficulty, that for p = + the operators in the 
HPT become nonhcrmitian at < N > > s, while for 
p = — , though this is avoided, Gaudin Hamiltonian (1) 
becomes nonhermitian to maintain its integrability. We 
resolve this problem by taking the spin limit sq — > oo in 
HPT for the radiation field at n = 0, yielding 

S+ = h\ S " = -b, Sl = -± e =-L^0. (2) 

We have considered here p = +1 for definiteness and 
kept terms up to order 0{\). For modeling N a two- level 
atoms we take spin sj — \ representation at all other 



Wl (0j) - 1, w 2 (0j) - k, w 3 (0j) - 0, w a (jk) 



This shows that in the given order now we have 
w-(jk) = 0, wa(j0) = 0, which simultaneously removes 
both the above obstacles and derives finally from (1) in 
the first nontrivial order O(-), the integrable matter- 
radiation models with mutually commuting Hamiltoni- 
ans [Hj, Hk] = 0: 

H j =w a oS+Hf r + HJ' r , (4) 

with j = 1, . . . , N a , where 

E h ° = n+(ba+ + ftto--) + n_(&t CT + + baj), (5) 

with Q± = 1 ± k , describes interaction between atoms 
and the radiation with explicit CRW terms and 

H r = E TT^rM^ + ^< + (6) 



points j — 1,2, ... , N a . The i coefficients in the expan- 
sion of the elliptic Gaudin Hamiltonians, though have the 
desired CRW terms, not yet yield the mutually commut- 
ing conserved set. The reason for this is twofold: the 
appearance of xyz type spin terms c^cr^ + cc. with co- 
efficients W-(jk) = ^(wi(jk) — W2(jk)) and the inhomo- 
geneity ws(j0) in the a 3 term, both of which we have to 
tackle before extracting integrable models from such an 
expansion. Obviously the first difficulty can be trivially 
resolved by setting k = 0, which yields W-(jk) = as 
considered in [16, 17]. However, since our aim is to keep 
the elliptic modulus k nontrivial, we have to take a dif- 
ferent rout. Our strategy would be to push the undesired 
w_ (jk) term out from the given order by suitably scaling 
the inhomogeneity parameters as z = K + ex , zj — 
exj, j = 1, . . . N a and taking the limit e — > 0, K be- 
ing the elliptic integral [19]. Observing that su{2) is 
isomorphic under the reflection of any axis of the basis 
vectors, we also redefine the coupling constants in (1) as 
w 2{jk) — > — W2(jk), Wi(jk) <-» w 3 (jk), and similarly for 
w a (0j) , to have more conventional notation. The rede- 
fined coupling constants using w a (0j) = —w a (j0) and the 
property of the elliptic functions, reduce in the needed 
order to 



*- -, a = 1,2,3, with k = (1 - k 2 )i. (3) 

e{ Xj -x k ) 

I 

models interatomic interactions, but without having any 
xyz type spin term. This crucial fact permits us to in- 
clude the atomic inversion term with arbitrary coefficient 
uj a in (4) , since ct| commutes with the whole Hamiltonian. 
The coupling constant f2_ for the interaction with CRW 
terms clearly vanishes at k = 0, recovering the earlier re- 
sults with RW approximation [16, 17]. Note that, we can 
construct a series of integrable Hamiltonians by different 
combinations of the commuting set Hj. For example, a 
generalization of the Tavis-Cummings model [8] without 
RW approximation can be constructed as 

3 3 

+ n_(&Vt + 6o7)], (7) 

without having explicit interactions between the atoms. 
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At k = 1, when Q + = 0_ = 1, (7) reduces further to 

H x = £> o0 f + (6 + 6t)(a+ + <77)), (8) 

which appears in the trapped ion model interacting with 
the center of mass motion, after taking the customary 
Lamb-Dicke limit. Notice that (8) is similar to the model 
studied in [12], if a bosonic number operator term is 
added to it. However, while the model in [12] becomes 
integrable only at the thermodynamic and the mean field 
limit of the atoms, our model (8) achieves this without 
going to such limits. Moreover, contrary to the pop- 
ular belief, that such a model with finite N a atoms is 
exactly solvable only under RW approximation and oth- 
erwise chaotic [13], we show it to be Bethe ansatz solvable 
in general, at least for w/ = 0. 

Since our model is derived from the elliptic Gaudin 
model through some limiting procedure, we can obtain 
its Bethe ansatz formulation also from the related result 
[15], by taking suitably the limits of the parameters in- 
volved. One important consequence of this inheritance 
is that, the excited Bethe eigenstates \M > , as in the 
elliptic Gaudin model, are no longer arbitrary, but con- 



strained by the total spin value M — J2 n s n . For our 
model therefore, we must have: M = \N a + s — > 1/e 2 
due to the limit s = 1/e 2 — > oo for the bosonic mode. 
This macroscopic excitation with high photon number is 
in conformity with the results of [12], indicating that we 
must be in the super-radiant phase with no possibility 
of phase transition, since due to uif = the critical cou- 
pling here would be ^/w/Wj = 0. Inspite of the fact that, 
the investigation of [12] is valid only for nontrivial u>f 
and its approach is totally different from ours, some ap- 
parent similarity between these results in the integrable 
case, perhaps may be explained by the point that of the 
thermodynamic limit adopted in [12] is mimicked by the 
high spin limit for the bosonic mode, considered in our 
model. 

For constructing exact eigenstates and eigenvalues, we 
introduce scaling also for the two sets of Bethe param- 
eters: Wb — eh,b — 1, ...^ and = el° a + Ka — 
1, . . . , s , similar to those for the inhomogencity parame- 
ters defined above. The exact eigenvalue for our matter- 
radiation models (4) is then derived from the limiting 
values of the Bethe ansatz result of [15] as 



where S = e 2 l° a . Possible nontrivial contribution in 
the given order may be obtained for S as S — 1°, when 
Z" = 1° are degenerate for all a = 1, ... , s , which wc 
consider here for defmiteness. From the Bethe equations 
for Z° we find a solution for these degenerate parameters 
as Z° — Xo, while for the rest of the Bethe parameters 
lb, b — 1, . . . , one needs to solve the equations 

5£ra-«- + t(fc^)- (10 » 

for a given set of arbitrary inhomogeneity parameters 
Xj,j = 1, . .., N a . For practical relevance we have ap- 
plied this Bethe ansatz result to a N a = 2 ion model 
with CRW terms described by (4) and shown its exact 
energy spectrum for j = 1 in Fig. 1. Similar result can 
be obtained easily for j = 2, by interchanging the param- 
eters x\ <-» x 2 - 

For a closer contact with physical systems we may in- 



clude a radiative excitation term ujfN, perturbatively, 
over our Bethe ansatz solvable models (4,7,8). For 
smaller values of field frequency ujf, when the RW ap- 
proximation worsens, the perturbative treatment would 
become more accurate. One can generalize the formu- 
lation of standard perturbation theory for applying it 
to integrable models with exact eigenvalues E n {l) like 
(9) and the corresponding normalized Bethe eigenstates 
with Bethe parameters I = (Zj, . . . , l® , l\, . . . , Im„ ). 
The result would however depend on the concrete model 
and might be difficult to extract in the explicit form. For 
considering the bosonic number operator term we are in- 
terested in, denoting N(l,m) =< ip(l)Ntp(rh) >, we can 
derive the first order perturbative correction to the exact 
eigenvalues as u)fN(l,l) and to the corresponding eigen- 
states as ujf - ?l} m ' , ^. ipim), where the sum is 

"l m J E n (l) — E n (m) v ' 

over the whole solution set of the Bethe parameters. 
Thus we have constructed and exactly solved multi- 
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important field excitation term perturbatively, over the 
exactly solvable models. 
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FIG. 1: Vacuum E0 and excitation energies E± showing 
asymmetric Rabi-type energy splitting, for different values of 
elliptic modulus k £ [0, 1]. The graphs depends crucially on 
the inhomogeneity parameters (chosen here as 2:0 = l,xi = 
1.1, xi = 3) and the atomic frequency (taken as u a = 1). 
However the qualitative behavior appears to remain the same: 
all the energies rapidly increasing with the increase of k, with 
a curious dip toward the end. 



atom matter-radiation models without rotating wave ap- 
proximation and with explicit interatomic interactions. 
We derive our intcgrable models and the corresponding 
exact Bethe ansatz result from the elliptic Gaudin model 
through a limiting procedure. We can include physically 
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